Charmed baryon S c (2800) as a ND hadronic molecule 
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The isotriplet E c (2800) baryon with possible quantum numbers J p = | or | is considered as 
a hadronic molecule composed of a nucleon and a D meson. We determine the strong two-body 



decay widths E c — > A c ty which are shown to be consistent with current data for the J p = | + and 
J p = | assignments. 
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A few years ago the Belle Collaboration [1[ observed an isotriplet of new baryon states with open charm S c (2800) 
decaying into A c 7r. This resonance was fit by a D-wave Breit-Wigner distribution based on the measured mass of 
the state and the Belle Collaboration tentatively assigned the quantum numbers J p = | . The same neutral state 
E° was possibly also observed in B decays by the BABAR Collaboration Q. Although the measured width of this 
resonance is consistent with the Belle value, the mass value is higher and somewhat inconsistent with the previous 
measurement. The BABAR Collaboration indicates that there is weak evidence that the excited they observe has 
| J = 1/2. In the following we will assume that both collaborations observe the same baryon resonance, although the 
■ present mass discrepancy and the final assignment of quantum numbers remain to be resolved. 

Earlier quark model predictions [3] for excited baryons containing one charmed quark lead to a possible identification 
of this state as a member of the J p = 3/2~ and 5/2~ doublet, where the light quark subsystem carries total angular 
momentum of 2 units. Further quark model studies on the mass spectrum of excited S states were also performed 
T — I , in [2-0] . These later works also tend to identify the S c (2800) as one of the nearly degenerate orbital excitations with 
L" ■ J p — l/2~,3/2~ or 5/2 _ . The strong decays of excited charmed baryons have been considered in the framework 
. £h ! °f heavy hadron chiral perturbation theory (HHChPT) Q and in quark models [H, [T3| • The computed A c 7r decay 
widths H of various J p assignments for the E c (2800) in the context of the 3 Pq model are found to be inconsistent 
with observation. In the chiral model of [l(| the A c 7t decay is consistent with a J p = 1/2 - assignment, while the 
HHChPT framework of Q uses the total width of E c (2800) as an input. Hence from the theoretical side a unique 
interpretation of this resonance in terms of an orbital excitation of the charmed three-quark system is presently not 
available. 

Alternatively, in Ref. [ll[ it was suggested that the S c (2800) is a so-called chiral excitation of open-charm baryons 
with J p = 3/2 _ - a charmed baryon resonance generated by S-wave coupled-channel dynamics of Goldstone bosons 
with the J p = 3/2 + sextet baryons with open charm. When including further pseudoscalar mesons like D mesons in 
the coupled-channel dynamics the updated results of [l2[ resulted in J p — 3/2~ state, where the width is much too 
large to justify the identification with the S c (2800). The dynamical generation of charmed baryon resonances in the 
context of a unitarized meson-baryon coupled channel model was also pursued in Refs. [HI, [l4 |. Both pseudoscalar 
and vector mesons are included in the S-wave coupled channel formalism as required by heavy quark symmetry. Now 
the S c (2800) is identified with a dynamically generated resonance in the J p = l/2~ channel with a dominant ND 
configuration. But the resulting partial A c 7r decay width is much too small to justify this identification. 
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Given the sparse experimental information concerning the excited E c spectrum and decays, but also the somewhat 
contradictory theoretical interpretations it is intruiging to note that the observed E c (2800) states are very close to the 
respective ND thresholds. For example, the measured mass difference Am = m s o( 2 8oo) ~ ™a+ = (515 ± 3^g) MeV 
by the Belle Collaboration [l[ corresponds to an absolute mass of ?7j£(>(2800) = (2802^) MeV [H| which should be 
compared to the nD° threshold value of about 2804 MeV. The closeness of these thresholds could imply that the ND 
components play a dominant role in the E c (2800) configurations, either by coupling of the excited three-quark state 
to the ND channels or in a hadronic molecule configuration. Although a full dynamical calculation was not performed 
yet concerning binding in the ND channel both for J = 1/2 and J — 3/2 (values for J as suggested by experiment) 
here we pursue a possible hadronic molecule interpretation of the E c (2800) - bound state of the nucleon and the 
charm D meson. Our aim is to work out in a hadronic framework the resulting E c — > A c ir decay widths for possible 
quantum number assignments of J p = l/2 ± and 3/2 ± which will be confronted with the experimental results. Note 
that while J p = l/2~ corresponds to an S-wave ND configuration, the options J p = 1/2+ or J p = 3/2 + represent 
a P-wave and J p = 3/2^ a relative D-wave in the ND system. Although slight binding in the ND seems less likely 
for higher par tial waves, especially for the D-wave, the possibility of such a weakly bound system is not excluded yet. 

InRefs. 161-1201 we developed the formalism for the study of recently observed exotic meson states (like D* (2317), 
D s i(2460), V(3872), V(3940), V(4140, •••) as hadronic molecules. The extension of our formalism to baryonic 
molecules has been done in Refs. [2l|,[22|]. A composite structure of these molecular states is defined by the compos- 
iteness condition Z = [23T - [25| (see also Refs. [HI]- [HI)- This condition implies that the renormalization constant of 
the hadron wave function is set equal to zero or that the hadron exists as a bound state of its constituents. The com- 
positeness condition was originally applied to the study of the deuteron as a bound state of proton and neutron [23| ■ 
Then it was extensively used in low-energy hadron phenomenology as the master equation for the treatment of mesons 
and baryons as bound states of light and heavy constituent quarks (see e.g. Refs. [U, HU). By constructing a phe- 
nomenological Lagrangian including the couplings of the bound state to its constituents and the constituents with 
other particles we calculated one-loop meson diagrams describing different decays of the molecular states (see details 
in [Hi-[22j]). 

In the present paper we proceed as follows. First, we discuss the basic notions of our approach. We consider a 
choice for the effective meson Lagrangian for the treatment of the E c baryons as ND bound states: E+ + = (pD + ), 
E+ = (pD° + nD + )/y/2, E° = (nD°), Second, we consider the two-body hadronic decays E c — > A c + ir. Finally, we 
present our numerical results. 

We consider the triplet (E+ + , E+, E°) as molecular states composed of nucleons and D mesons as: 

|E++) = \pD+), 

|E+) = -L\ p D° + nD + ), (1) 
|E°) = \nD°). 

Our approach is based on an effective interaction Lagrangian describing the couplings of the E c to its constituents. 
The molecular structure of the E c baryon with quantum numbers J p = | is described by the Lagrangian 

£ Sc (a:) = <? Ec S;(z)J Sc (z) + H.c, Je c {x) = D{x) t T J d 4 y $(y 2 ) N(x + y) (2) 

while for the choice J p = | the Lagrangian contains a derivative TV D coupling 

C^(x) = 9sc ££(z)J Ec>M (z) + H.c, J E » = D(x)rr f d 4 y^(y 2 )d^N(x + y) (3) 



where g s is the coupling constant of the isotriplet S c to the N = (p, n) T and D = (D°, D + ) T constituents. Here T is 

1+3- 
2 ' 2 



the corresponding Dirac matrix related to the spin-parity of the E c . In particular, for J p = i , | we have T = 7' 5 



while for J p = | , | + the Dirac structure r = / should be inserted in the E c current. 

We propose a picture for the E c in analogy to heavy quark-light antiquark mesons, i.e. the heavy D meson is 
located at the center of mass of the E c , while the light nucleon surrounds the D. We describe the distribution of the 
nucleon around the D meson by the correlation function ^(j/ 2 ) depending on the relative Jacobi coordinate y. A basic 
requirement for the choice of an explicit form of the correlation function $(y 2 ) is that its Fourier transform vanishes 
sufficiently fast in the ultraviolet region of Euclidean space to render the Feynman diagrams ultraviolet finite. We 
adopt a Gaussian form for the correlation function. The Fourier transform of this function is given by 

$(p|/A 2 ) = exp(-pi;/A 2 ) , (4) 
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FIG. 1: Diagram describing the E c mass operator. 




FIG. 2: Diagrams contributing to the E c — > A c 7r decay. 



where pe is the Euclidean Jacobi momentum. Here, A ~ m^v ~ 1 GeV is a size parameter, characterizing the 
distribution of the nucleon in the E c baryon, which is of order of the nucleon mass or 1 GeV. In the numerical analysis 
we therefore fix the mean value to A = 1 GeV. 

The coupling constant g Sc is determined by the compositeness condition [l6|, I23U25] ] . It implies that the renormal- 
ization constant of the hadron wave function is set equal to zero with: 

Z Sc = 1-E'sJmsJ =0. (5) 

Here, £^("1^) = 9^ (m^J is the derivative of the mass operator for J = |. For J = | the same relation holds 
but now Ss c (msJ should be identified with the scalar function proportional to the Minkowski tensor g^ v in the full 
mass operator ££". Note, that for J = | the other possible Lorentz structures in the E c mass operator vanish due to 
the Rarita-Schwinger conditions. The mass operator of the E c baryon is described by the diagram of Fig.l. To clarify 
the physical meaning of the compositeness condition, wc first want to remind the reader that the renormalization 

1/2 

constant Z E can also be interpreted as the matrix element between the physical and the corresponding bare state 
— an elementary structureless field. For Zs c = it then follows that the physical state does not contain the bare one 
and hence is described as a bound state. As a result of the interaction of the E c baryon with its constituents N and 
D, the E c baryon is dressed, i.e. its mass and its wave function have to be renormalized. Note, in the present paper 
we only consider the contribution of a possible molecular (ND) component to the structure of the E c . An inclusion 
of a three-quark component is possible, but goes beyond the scope of the present paper. 

In Table I we display the results for the coupling g + of the single charged E+ state for different spin-parity 

assignments and for a variation of the size parameter A in the region of 0.75 - 1.25 GeV. Note, that an increase of the 
A value leads to an enhancement of the couplings g + . A final value for the cutoff model parameter A can ultimately 

only be fixed when more decay data on E c are available. 

The one-loop hadron diagrams contributing to the E c — > A c + it transition are displayed in Fig. 2(a) and 2(b). To 
evaluate these strong E c — > A c + n decays we further need an effective Lagrangian including the couplings of the E c 
constituents to A c and ir. The irNN and D*Dtt couplings are constrained by data or by low-energy theorems. In 
particular, the nNN coupling at leading order of the chiral expansion is expressed through the nucleon mass m^r, the 
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Table I. Coupling constants g s + for different J p assignments. 
Error reflects variation in A from 0.75 to 1.25 GeV. 



T P 1 + 

J ~ 2 


J ~ 2 


tP 3 + 

J ~2 


J ~2 


7 ± 1.9 


0.6 ± 0.2 


4.2 ± 1.4 GeV" 1 


35.3 ± 1.8 GeV -1 



pion decay constant = 92.4 MeV and the nucleon axial charge §a = 1.2695 as: 

£~kNN = 9-kNnN vysn r N , g-^NN = ~et- 9 a ■ (6) 

For the D*Dir coupling we take the central value of gD*D-n — 17.9 extracted from the measured D* + — > D° + tt + 
decay width [26| 

T(D*+ -> £)°7r + ) = 9 ^*^ P* 3 , (7) 
247rm^, + 

where P* is the three-momentum of 7r + in the D* + rest frame. Then the interaction D*Dtt Lagrangian reads: 

C D , Dv = ^^D*Ud^TD + H.c. (8) 
v2 

The couplings DNA C and D*NA C are estimated by matching the flavor SU(4) effective Lagrangian to the SU(3) 
version, both describing the couplings of pseudoscalar and vector mesons to two baryons (see details in the Appendix): 

9DNA C — —g-irNN , gD»NA c — 9pNN ■ (,»J 

In the evaluation of the diagrams in Figs.l and 2 we use the standard free propagators for the intermediate particles: 
iS N (x -y) = (0\TN(x)N(y)\0) = [ e^-"^*), S N (k) = 1— — (10) 



for nucleons, 



d A k 1 
iS D (x -y) = (0\TD(x)D\y)\0) = / — — e~ ik ^S D (k) , S D (k) = — -_- (11) 



(2tt) 4 « " y " UK ' m 2 D -k 2 -i 



for pseudoscalar D mesons and 



iS%(x-y) = (0\TD*»(x)D**(v)\0) = / ±^ S^,{k) , S%(k) = 9 2 + \ 2 ' D ' 

x ' J (2ir)^i vn D , — K — id 



(12) 

j y&n j ~l iiljj* — — it 

for vector D* mesons. 

The decay widths of the strong two-body transitions S c — > A c + it are then calculated according to the following 
expressions for the different spin-parity assignments S C [J P ]: 

r(E c [l/2+]^A c + 7 r) = -^^A 1 / 2 (m| c ,mi,m2 ) ( (mSc _ mA j2_ TO 2^ (13a) 
r(E c [l/2-]^A c + 7 r) = ^^xy 2 (ml c ,ml c ,ml)((m Sc +rn Ac ) 2 -mi), (13b) 
r(E c [3/2+]^A c + 7 r) = A 3 / 2 (m| e ,m 2 e ,m 2 )(W -m A J 2 -ml) , (13c) 

2 

r(E c [3/2-]^A c + 7r) = g' A °* A 3 / 2 (m|. , m\ c , ml) ((m s „ + m Ae ) 2 - m 2 ) , (13d) 
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Table II. Effective couplings <fe c A C 7r, /is c A c7r , /s c A c7 r and gv c A c w 
Error reflects variation in A from 0.75 to 1.25 GeV. 



Mode 


d£ c A c 7r 


^S c A c 7r 


/e c A c tt 




£++ -> A+7T + 


- 8.15 ± 2.72 


1.63 ± 0.54 


1.95 ± 0.97 GeV" 1 


- 3.35 ± 1.61 GeV" 1 


E+ -> A+tt 


- 7.78 ± 2.60 


1.48 ± 0.47 


1.90 ± 0.95 GeV" 1 


- 3.24 ± 1.59 GeV -1 


E° -»■ A+TT" 


- 7.52 ± 2.54 


1.43 ± 0.45 


1.87 ± 0.94 GeV" 1 


- 3.16 ± 1.58 GeV" 1 



Table III. E c — > A c tt decay widths (in MeV) for different spin-parity assignments of the E c . 
Error reflects variation in A from 0.75 to 1.25 GeV. Results for preferred value of A = 1 GeV are given in brackets. 



Mode 


J " 2 




tP 3 + 

J — 2 


jP _ 3" 
■ J — 2 


£++ A+7T+ 


41.1 ± 24.7 (37.0) 


173.6 ± 103.6 (156.4) 


0.176 ± 0.140 (0.141) 


54.5 ± 42.6 (44.3) 


E+ A+TT 


37.6 ± 22.6 (33.9) 


142.3 ± 82.1 (129.3) 


0.171 ± 0.137 (0.137) 


51.9 ± 41.0 (41.8) 


E° -> A+7T- 


35.1 ± 21.3 (31.5) 


132.3 ± 75.8 (120.4) 


0.164 ± 0.132 (0.131) 


49.5 ± 39.6 (39.6) 



where X(x, y, z) = x 2 + y 2 + z 2 — 2xy — 2yz — 2xz is the Kallen function; ms e , m\ a and are the masses of S c , A c 
baryons and the pion. We also introduce the effective coupling constants d£ c A c7r , ft.s c A c ir /s A C 7r and 5s c a cT defining 
the interaction of the S c of definite spin-parity with A c and n as a result of the processes in Fig. 2 with: 



£s c (i/2+)A c7 r = ds c A c 7r A c 75 7r S c + H.c. , (14a) 

£ Sc (l/2-)A c7 r = ^S c A c7 rA c 7rS c + H.C. , (14b) 

^S c (3/2+)A c ,r = /EoAcTrAT^TrS^ + H.C., (14c) 

^e c (3/2-)a c tt = .gs c A c 7r A c 7 5 a M 7r + H.c. . (14d) 



In Table II we present our results for these effective couplings <fe c A e7r , ft.s c A c 7r, /e c a c 7t and <7s c a c 7t including a variation 
of the cutoff parameter from 0.75 to 1.25 GeV. 

Our final numerical results for the decay widths are summarized in Table III. For the E c masses we use the 
measured values [l| of the Belle Collaboration. The predictions for the decay widths differ sizably depending on the 
J p assignment for the S c (2800). These predictions are to be compared to the measured total widths of the S c (2800) 
baryons [lj,ll5j with: 

r(S++) = 75li!±u MeV, 

r(S+)=62lg+52 MeV , (15) 
r(S°) = 61±\ s +il MeV. 

Since the observed A c 7r decay modes of the S c (2800) states are assumed to be dominant present results favor, at least 
in the context of the ND molecule interpretation, either the J p = 1/2+ or the J p = 3/2~ assignment. Note, the 
J p = 3/2 _ assignment was originally assumed by the Belle Collaboration [l|. The alternative scenario for E c with 
J p = 3/2+ is clearly excluded by the predictions of Table III. For small values of the dimensional parameter A the 
scenario for S c with J p = l/2~ becomes compatible with data. 

In conclusion, we estimated the strong A c tt decays of the S c (2800) baryon for different spin-parity assignments 
assuming a dominant molecular ND structure of this state. Judging from the decay widths of the order of 40 MeV 
we find that the original scenario where the S c has spin-parity J p = | and the choice J p = ^ + are consistent 
with current data. The option J p = | + leads to strongly suppressed partial decay widths of the order of a hundred 
keV, while J p = \ leads to enhanced partial decay widths and only becomes compatible with data for relatively 
small values of the dimensional parameter A. Although weak binding in the ND system for J p = 1/2+ and J p — 
3/2 ~ remains to be studied, present evaluation of the A c 7r decay widths point to a possibly sizable role of the ND 
configuration in the E c (2800). 
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Appendix A: Matching of the phenomenological SU(3) and SU(4) PBB and VBB interaction Lagrangians 

First we consider the PBB interaction that is the coupling of a pseudoscalar (P) meson to two baryons (BB). In 
flavor SU(3) the couplings are generated from the 0(p) term of chiral perturbation theory (ChPT) [27] describing the 
coupling of baryon fields with the chiral fields: 

£ S p U b 3 b = --^J^ mB + m B )tr[B tl 5 {PB}) - ^j={m B + m s )tr(Bi^[PB]) . (Al) 

Here Fp = F n = 92.4 MeV is the leptonic decay constant; D and F are the baryon axial coupling constants (we 
restrict to the SU(3) symmetric limit, where D = 3F/2 = 3(^/5 with = 1.2695 being the nucleon axial charge) 
the symbols tr, {. . .} and [. . .] denote the trace over flavor matrices, anticommutator and commutator, respectively. 
We replace the pseudovector coupling by the pseudoscalar one considering on-mass-shell baryons. The SU(3) baryon 
B and pseudoscalar meson P matrices read as: 



B 



( ^°/V2 + A/VQ £+ p \ 

£- -£°/V2 + A/V6 n 

5- S° -2A/V6 J 



(A2) 



P 



7T+ K+ \ 

-ir /V2 + v /V6 K° . (A3) 
K° -2??/V6 / 

The SU(4) PBB Lagrangian is given by [26 

C S P % = 9l B kmn i lb P l k B lmn + g 2 B kmn i lb P l k B lnm , (A4) 




where the indices l,m,n of the tensor Bi mn run from 1 to 4, representing the 20-plet of baryons (see details in 

•Dl 

k 



Refs. 28]); P l k is the matrix representing the 15-plet of pseudoscalar fields. The baryon tensor satisfies the conditions 



(A5) 

The full list of physical states in terms of SU(4) tensors is given in Ref. (28|. Here we only display a few of them: 

V = -B112 = — 2.B121 = — 2£?2ii , n — — B221 = 2i?2i2 = 2£?i22 , 

E+ + = Bn4 = — 2_Bi4i = — 2S 4 ii , S° = — B 2 24 = 2i?242 = 2B422 , (A6) 
7T+ = Pi , Tr - = P\ , D° = Pi , D* + = Vi , D*° = Vl . 

Evaluating the irNN couplings in both versions we fix the SU(4) couplings g\ and g 2 as (in the SU(3) Lagrangian we 
restrict to the mass degenerate case ms = nig = m p — 938.27 MeV): 

5 D-F gi + g2 . 

9-knn = gi - -^92 , g-KNN D + F = — 4 • ( A7 ) 

Considering the SU(3) symmetric ratio of F and D couplings F/D = 2/3 we get 

4 r 

9i = , .92 = --V2g 7lNN . (A8) 
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Finally the gDNA c coupling is fixed as 

9dnk c = —g-xNN ■ (A9) 

In complete analogy we fix the vector meson VBB couplings. The SU(3) VBB Lagrangian can be expressed in terms 
of the pNN coupling constant as: 

£ s v u i B = + ^tr(VsK, (A10) 

where 

/p°/V2 + ^/V2 p+ _ K*+\ 

V=\ p- -p°/V2 + io/V2 K*° . (All) 

\ K*- K*° -<j) J 

The SU(4) VBB Lagrangian is given by [28) : 

= h x B kmn ^Vl k B lmn + h 2 B kmn ^V^ k B lnm . (A12) 

Evaluating the pNN couplings in both versions we fix the SU(4) couplings hi and h% as: 



Finally, the D*NA C coupling is fixed as 



hi = 2h 2 = ^~^9 p nn ■ (A13) 



V3 

9D*NK C = 7^9pNN , (A14) 



where for the g P NN coupling we take the SU(3) prediction of 

g P NN = 6 . (A15) 
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